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Abstract. Earlier the quasi-isothermal Stackel’s model of mass distribution 
in stellar systems was considered. The expression for spatial density was found. 
In this work an application of this model to our Galaxy is studied. The model 
rotation curve was htted to data on kinematics of neutral hydrogen. Estimations 
of structural and scale parameters of the model and equidensities for our Galaxy 
are presented. 
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1. INTRODUGTION 

There are various methods of modelling the gravitational potential of our 
Galaxy and other stellar systems. One of them is based on assumption that the 
potential is of Stackel’s form. After pioneering works by Kuzmin (1952, 1956) the 
method was further developed by de Zeeuw et al. (1986), Dejonghe & de Zeeuw 
(1988) and others. In this paper we study possibilities of using the quasi-isothermal 
potential suggested by Kuzmin et al. (1986) for constructing Stackel’s model of 
our Galaxy. 

It is known that the third quadratic in velocities integral of motion 

h = [Rvz - ZUR^ + z^vl P zl{vl - 24>*), (1) 


exists for Stackel’s potentials. Here R, z are the cylindrical coordinates; vr, vg, 
Vz are the projections of spatial velocity; zq is a scale parameter of dimension of 
length; a function ^*{R,z) must satisfy the equations 
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where 4)(R, z) is an axisymmetric potential. The assumption on the existence I 3 
allows to explain the observed triaxiality of velocity ellipsoid. 
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In this work we suppose that the potential in the equatorial plane has the 
following form: 


= $oln 


^ w{R) 


( 3 ) 


where j3 G [0, +oo) is a structural parameter of the model, 


w^{R) = 1 + k^R'^, 


( 4 ) 


<i>o and K are scale parameters. 

This potential was proposed by Kuzmin et al. (1986) for spherical systems and 
was called quasi-isothermal. The potentials by Schuster-Plummer and Jaffe are 
limiting cases of this one (when (3^0 and /3 —>■ oo respectively). 

Stackers models of mass distribution with the quasi-isothermal potential were 
constructed and an analytical (though cumbersome) expression for a spatial den¬ 
sity was found (Gromov 2012, 2013, 2014). Stackel’s models with the Schuster- 
Plummer and Jaffe potentials were also studied (Gromov 2013, 2014). 

2. DATA 

In this study we use data on the rotation of the neutral hydrogen: hve inde¬ 
pendent data sets from the whole 21-cm line prohle and one set from the 21-cm 
tangent points (see details in Nikiforov & Petrovskaya 1994). The whole profile 
data contain Gamm’s function values = Ro{ui — wlsr) in relation to x = R/Rq, 
where R is the distance to the galactic axis, ui and wlsr are the angular velocities 
of HI at the distance R and of the Local Standard of Rest respectively, Rq is the 
distance of the Sun to the Galactic center. The tangent points’ data also define in 
fact the dependence of 12 on x. The total number of HI data points is 239. 

It should be pointed out that analyses of the whole 21-cm profiles give at any 
radius the velocity of rotation averaged over essentially all galactocentric circle 
of this radius. Therefore the resultant HI rotation curve is well smoothed over 
essentially all galactocentric longitudes. 

A value of the linear velocity of rotation, Vi, for an HI data point can 

be calculated by the formula 


Xi = -I- xlsr] a:i, (5) 

where ulsr is the linear velocity of the LSR. Here we adopted a value for ulsr = 
220 km/s as a combination of constants ulsr = ^lsR) where Rq = 8.0 kpc (e.g., 
Reid 1993; Nikiforov 2004; Nikiforov & Smirnova 2013), wlsr = 27.5 kms“^ kpc“^ 
is an intermediate value between recent determinations (e.g.. Feast & Whitelock 
1997; Zabolotskikh et al. 2002). Since as a rule the number of points in a HI data 
set can be arbitrary set by authors of study, the weights of data points, pi, were 
adopted to be proportional to the length of interval [xminj a^max] covered by this 
set (see Table 1 in Nikiforov & Petrovskaya 1994). 

Notice that the form of the rotation curve based on the HI data is independent 
of calibrations of distance scales and is specified by only the parameter xlsr, 
according to Eq. ®. Choosing a value for Rq affects the value of xlsRi with our 
parametrization of the HI rotation curve, and scales the coefficient k in Eq. 
in any case. Other parameters of the model potential ([3]) do not directly depend 
on an adopted value of Rq. 
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Fig. 1. Comparison of the model circular velocity curve with data on the rotation of 
the neutral hydrogen. 


3. ESTIMATION OF $o, P, k 

To construct the Stackel’s model potential ([3]), (|4]) for our Galaxy it is necessary 
to estimate parameters /?, k. We shall fit the model circular velocities, Vc, to 
observational data on the rotation of the Galaxy. An expression for the circular 
velocity 

v!{R) = -R^iR,0) ( 6 ) 

in the case of the quasi-isothermal potential has the following form: 


v^,{R) = d>o- 

{1 + k^R^) 


flT , ^ 

V VTTk^ 


( 7 ) 


where q 


P 

/3 + 1 ' 


We shall use g as a structural parameter instead of /3 for 


ge [0,1]. 

The model parameters were estimated by the ordinary least-squares fitting. 
We minimized the statistics 


239 

= '^Pi [vc{Ri) - Vi]'^ , 
2 = 1 


( 8 ) 
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where Vc{Ri) is the model value of circular velocity at Ri calculated by Eq. ([7]), Vi 
is the “observed” value of circular velocity calculated by Eq. m, Pi is the weight. 

We found, that minimum is achieved, when $o = 258.1 ± 1.5 km^s“^, 
K = (0.32 ± 0.01) {Ro/8 kpc)”^ kpc“^, q = IIq.oos • With these values, the quasi- 
isothermal model gives the best approximation for observational data. 

A comparison of the model curve of circular velocity with observational data is 
shown on Fig. 1, where the solid curve is the model velocity curve vdR) and points 
are the HI data. The mean error of weight unit for this solution is ct = 2.98 km s“^. 
Although a polynomial model of rotation curve can reveal a finer structure of the 
HI rotation law with a = 2.10 kms“^ for the same data (Nikiforov 2000), Fig. 1 
demonstrates that the quasi-isothermal potential in the equatorial plane represents 
the main trend in the rotation low and can be used for an approximation of the 
rotation curve for our Galaxy. 

It is worth noting that the form of rotation curve for the quasi-isothermal po¬ 
tential (Fig. 1) is similar to the form of the “universal” rotation curve for spiral 
galaxies (Persic et al. 1996). Recent results of kinematic modelling the Galactic 
subsystem of high-mass star forming regions (Reid et al. 2014) demonstrate that 
the universal rotation curve of Persic et al. (1996) is well suited to the representa¬ 
tion of the rotation curve of the Galaxy. 


4. ESTIMATION OF zq 

The expression for density of Stackel’s models includes the parameter zq. But 
it is impossible to find the latter from the rotation curve. In principle, zq can 
be estimated from the run of total spatial density, but presently we cannot do 
it for any one-component Stackel’s model. Indeed, we do not know the density 
run for the dark matter from observations. Hence, to estimate zq we are obliged 
to suppose that in the Solar neighborhood the potential of our model is close to 
potentials of other models constructed with using not only the rotation curve but 
other observational data as well. The model by Gardner et al. (2011) is the last 
of such models. Its potential is 


$ = -f + ^3 , 


(9) 


where 


<^>h = Iv^HR^ + z^ + RI), 

GMci GMc2 

-f -f Rq_^ R? + z'^ + Rq^ 




-GMd, 


1 ^R^ + {ad, + ’ 





where G is the gravitational constant, the values of parameters Vj^, Ri, 
Rci, Rc^i Aldi, adi, b, Mg„, bg are given by Gardner et al. (2011). 
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R, kpc 

Fig. 2. The function |2o(^)| for the model of Gardner et al. (2011). 


For any potential $(i?, z) it is possible to define a function 


z2(i?) 


{,d^R,z) , ^fd^^R.z) ^d^^R,z)\ 

+ dR^ dz^ ) 

d^<S>{R, z) 
dz'^dR 


( 10 ) 


(Kuzmin 1952, Ossipkov 1975). It follows from the condition of existence of the 
third quadratic integral (CD that Zq{R) = const for Stackel’s models. Fig. 2 shows 
|zo(i?)| for the model of Gardner et al. (2011). We see that this function is almost 
constant for R from 1.5 kpc to 9 kpc. The latter means that for such potential the 
expression CD can be used as an approximate integral of motion for stars moving 
inside this zone (with not very large l^zD- Functions zq{R) were constructed 
earlier by Einasto & Riimmel (1970) and Ossipkov (1975) for Einasto’s models of 
our Galaxy and M 31. These authors discussed the negativity of Zq{R) for large 
R that takes place also for the potential 

We conclude that it is possible to set zq = 5.4 kpc (this is a value of |zo| at 
i?o = 8 kpc for the model of Gardner et al. 2011). Earlier Ossipkov (1975) found 
zq = 7 kpc, and according to Kuzmin (1956) zq = 3.6 kpc. Our value of zq is close 
to zo = 4.8 kpc obtained by Malasidze (1973). 
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Fig. 3. Equidensities for obtained values of parameters. 

5. MODEL OF THE GALAXY WITH THE QUASI-ISOTHERMAL 
POTENTIAL 

For constructing the quasi-isothermal model of Galaxy the obtained values of 
parameters were substituted in an expression for spatial density (Gromov 2012, 
2013, 2014). As a result, the equidensities for central layers of Galaxy were ob¬ 
tained. They are shown on Fig. 3. As earlier (Gromov 2012, 2013, 2014), po is a 
density in the center of the model, when k = 1 kpc“^, /3 = 1, $0 = 1 km^ s“^. 

These equidensities are similar to equidensities for Jaffe’s model, which is the 
limiting case (/3 —>■ 00 ) of quasi-isothermal model 

4>(i?,0) = $oln(l-bA-iR-i) , (11) 

where A is a scale parameter. 

Fig. 3 shows a drastic drop in density at equidensities which dimensions are 
close to ones of the Galactic triaxial bulge/bar—(3.5 : 1.4 : 1.0) kpc (see Gardner & 
Flynn 2010). It is possible that this resemblance is accidental. However, the direct 
inclusion of bulge component in the dynamic model with the use of the same HI 
data (Nikiforov 2001) gives a similar result: a severely limited from above estimate 
of the cut-off radius of bulge Ub = 2 .II 2 U kpc. 

It is evident that the total mass of the quasi-isothermal model is 

M = ,12) 

The minimum value of M for the obtained values of parameters is 6 x 10^^ Mq . The 
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maximum value of the mass for the same parameter values is unlimited, because 
P can take infinity large values. 

6. CONCLUSIONS 

We estimated parameters of the quasi-isothermal model of mass distribution 
for our Galaxy using data on the rotation of the neutral hydrogen and obtained 
equidensities for this set of parameters. 

The comparison of the model with observational data shows the relevance of 
the model. So, the quasi-isothermal model can be used for approximation of the 
gravitational potential of our Galaxy. It is possible to apply the quasi-isothermal 
model to external galaxies with almost flat rotation curve. 

In the future a multi-component Stackel’s model will be constructed for the 
Galaxy to gain a better approximation to the data. 
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